Electronic structure of the two-dimensional Heisenberg antiferromagnet VOCl: a 

multi-orbital Mott insulator 
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We have studied the electronic structure of the two-dimensional Heisenberg antiferromagnet VOCl 
using photoemission spectroscopy and density functional theory including local Coulomb repulsion. 
From calculated exchange integrals and the observed energy dispersions we argue that the degree 
of one-dimensionality regarding both the magnetic and electronic properties is noticeably reduced 
compared to the isostructural compounds TiOCl and TiOBr. Also, our analysis provides conclusive 
justification to classify VOCl as a multi-orbital Mott insulator. In contrast to the titanium based 
compounds density functional theory here gives a better description of the electronic structure. 
However, a quantitative account of the low-energy features and detailed line shapes calls for further 
investigations including dynamical and spatial correlations. 
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PACS numbers: 71.20.-b,71.27.+a,71.30.+h,79.60.-i 



I. INTRODUCTION 

Transition metal compounds arc characterized by a 
plethora of many-body phenomena which emerge from 
a complex interplay of charge, spin, orbital, and lat- 
tice degrees of freedom. More recently, strong focus has 
been placed on the role of multi-orbital physics, which is 
governed by a competition between intra-atomic Hund's 
rule exchange coupling and crystal field splitting in the d 
shellfii^i^ The simplest possible case where multi-orbital 
effects may occur is in a (P configuration, as encoun- 
tered, e.g., in trivalent vanadium compounds. A promi- 
nent example is V2O3, whose spin nature, whether be- 
ing S' = 1, as naively expected for a d? system and 
found spectroscopically^^ or whether V-V pair dimeriza- 
tion leads to an effective S* = 1/2 behavior has been 
lengthly debated and recently settled as a 5 = 1 system 
(for details, see Ref. H). 

Here we study another prototypical V^"*" compound, 
namely VOCl. It belongs to the class of low- 
dimensional transition metal oxyhalidcs of the form 
MOX (M=Ti,V,Cr,Fe; X=Cl,Br,I). Amongst these, es- 
pecially the 3c?^ systems TiOCl and TiOBr have seen con- 
siderable interest recently^ii^iiSiiSiiii— as they also show 
dimerization leading to non-conventional spin-Peierls be- 
havior. Like its Ti counterparts, VOCl has a layered 
crystal structure which provides a topologically frus- 
trated lattice, thus being a good candidate for unusual 
ground states and phase transitions. Such tendencies 
are enhanced due to competing one- and two-dimensional 
interactionsj ^"^'^^i^^'^^i^^ Furthermore, it has been shown 
that TiOX can be electron-doped by alkali metal in- 
tercalation into the van der Waals gaps of the crystal 
structurei^ up to a nominal d^~^^ configuration. In some 
sense, the electronic structure of VOCl can be viewed as 




FIG. 1: (Color online) Left: Three-dimensional view of the 
VOCl crystal structure along the 6 axis. Two bilayers sep- 
arated by a van-der- Waals gap are shown. Top right: local 
octahedron around a V ion with the corresponding crystal 
axes. Lower right: triangular V sublattice projected onto the 
ab plane and viewed along c. Also shown are the occupied 
d^2_y2 (left) and d^z (right) orbitals to illustrate the hopping 
paths provided through overlapping electron clouds. In real- 
ity, neighboring V ions along a are shifted along c (see left 
structure). 



the X = 1 limit of such a doping series and thus provides 
an interesting reference point. 

Electronically, VOCl is found to be an insulator. Our 
photoemission spectroscopy (PES) data as well as the 
density functional theory (DFT) calculations confirm the 
importance of correlations in this system, and our anal- 
ysis shows that the suppression of metallicity is driven 
by strong correlations. This implies that VOCl is a 
multi-orbital Mott-Hubbard system, in line with results 
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FIG. 2: (Color online) Magnetic susceptibility of VOCl mea- 
sured with an external field of 0.1 T applied along the different 
crystallographic axes over an extended temperature range. At 
TVeei = 79 K a transition into an antiferromagnetic phase oc- 
curs, accompanied by an anisotropic behavior, evidenced by 
the sharp drop if the external field is applied along a. The 
inset shows specific heat measurements which display a peak 
anomaly at the same temperature. 

from TiOCl, which displays multi-orbital physics upon 
electron dopingi ^^i^°i^^ . On the other hand, the effects 
of dynamical and spatial fluctuations, which should be 
prominently observable in angle-resolved photoemission 
(ARPES), turn out to be less important than, e.g., in 
TiOCl, though not ncghgiblej^i^ 

This paper is structured as follows: in Sec.[II]we give an 
overview of the structural, thermodynamic, and magnetic 
properties of VOCl, in Sec. IfffI both experimental and 
computational methods are described. Section [IV] con- 
tains our results and a comparative discussion of the elec- 
tronic structure of VOCl at room temperature and above, 
determined using (AR)PES as well as DFT calculations 
within the generalized gradient approximation including 
local correlations at the mean field level (GGA-I-U). Fi- 
nally, in Sec.|V]we summarize our results. 



II. STRUCTURAL, THERMODYNAMIC, AND 
MAGNETIC PROPERTIES 

VOCl crystallizes in an orthorhombic structure of 
buckled V-0 bilaycrs sandwiched by CI layers as shown 
in Fig. [1] (space group Pmmn; lattice parameters a = 
3.78A, b = 3.30A, c = 7.9lA),^ These building blocks 
are only weakly coupled through van-der-Waals forces 
along the crystallographic c direction. Each V ion lies in 
the center of a strongly distorted octahedron of two CI 
and four O ions. Thus, the degeneracy within the low- 
lying t2g triplet and the high-lying Cg doublet is com- 
pletely hftcd, as is also the case for TiOCUiiii^ It will 



be shown below that the electrons occupy the lowest- 
lying orbitals of d^2_y2 and dxz character in the chosen 
reference frame (x = h,y = c,z = a), which leads to 
a weaker electronic anisotropy within the ab plane com- 
pared to TiOCl. This enhanced two-dimensionality is 
also reflected in the magnetic susceptibility displayed in 
Fig. O While showing an isotropic behavior for the high- 
temperature phase, a sharp drop with the external field 
applied along the crystallographic a axis is observed at 
around T^r^e; = 79 K, while along the b and c axis a slight 
upturn occurs. This behavior, in conjunction with the 
anomaly in the specific heat shown in the inset of Fig. [21 
indicates a transition to an antiferromagnetically ordered 
phase, which has been described previously by a two- 
dimensional S = 1 Heisenberg model on a square lattice 
with antiferromagnetic nearest-neighbor interactioni^ In 
this low-T Neel phase the crystal symmetry is reduced 
from orthorhombic to c axis monoclinicj ^^'^^ 



III. TECHNICAL DETAILS 

Single crystals were grown by chemical vapor 
transport Crystals of typical dimensions (3 x 2 x 
0.1) mm"^ were selected for the PES measurements. The 
samples were characterized by x-ray diffraction and mag- 
netic susceptibility measurements (see Fig. [5]). For the 
(AR)PES experiments we used a Specs PHOIBOS 100 
as well as an Omicron EA125 analyzer, monochromati- 
zed Al Kq radiation for x-ray photoemission spectroscopy 
(XPS; /ii/=1486.6cV) and Hel„ photons {hv^2l.2cY) 
from a gas discharge lamp for valence band spectroscopy 
(UPS, ARPES). The total energy resolution amounted to 
w700meV and 70meV for XPS and UPS, respectively. 
In ARPES, using the EA125, the angular resolution was 
about 1°. Clean surfaces were obtained by in situ cleav- 
age using Scotch tape. Their cleanliness and atomic long- 
range order were evidenced using XPS and low-energy 
electron diffraction (LEED), respectively. Pressure dur- 
ing measurements was in the low 10"^" mbar regime. Due 
to the strongly insulating nature of VOCl all measure- 
ments have been performed at temperatures of 360 K or 
above in order to minimize charging effects via thermally 
activated conductivityi^ It was found from systematic 
temperature variations that for T > 360 K the maximum 
of the V 3d peak saturates at ~2.2eV below the chemi- 
cal potential fj,exp as determined from the Fermi edge of 
a silver foil. All measured spectra have been aligned ac- 
cordingly, and the calculated spectra have been shifted 
such that their first order moments correspond to those 
of the experimental ones. 

The DFT calculations where performed with the 
Full Potential Linearized Augmented Plane Wave basis 
(FPLAPW) as implemented in the WIEN2k code^. We 
considered a k mesh of (15 x 17 x 7) in the irreducible 
Brillouin zone and a RKmax=7. For the GGA-I-U calcu- 
lations we considered the implementation of Dudarcv et 
alM- in WIEN2k where a unique parameter Ugf f = U ~ J 
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FIG. 3: (Color online) Comparison of calculated partial densi- 
ties of states (pDOS) for Ueff = 1.3 eV with angle-integrated 
PE spectra (no background subtraction). The former have 
been weighted with the respective photo-ionization cross- 
sections at the corresponding excitation energies (top: UPS, 
Helc,; bottom: XPS, AlKc). Also shown are the sums of 
the weighted pDOS (solid red curves) and their convolutions 
with a Gaussian in order to mimic the experimental energy 
resolution (grey-shaded areas). 



is needed as an input. This leaves a certain flexibility for 
the choices of U and J values provided their difference is 
Uef f ■ NMTO downfolding calculations^ were also per- 
formed in order to obtain the important hopping param- 
eters among the V atoms. 



IV. RESULTS AND DISCUSSION 

Figure [3] presents the angle-integrated UPS and XPS 
valence band spectra of VOCl together with the density 
of states (DOS), calculated for the room temperature 
structure within GGA-hU for Ueff = U - J = 1.3 eV 
where U is the onsite Coulomb energy and J the exchange 
integral. Also shown are the partial DOS (pDOS) for CI 
3p, O 2p, and V 3d which have been multiplied by cal- 
culated photo-ionization cross-sections^ for UPS (upper 
panel) and XPS (lower panel) excitation energies, respec- 
tively. In order to mimic the experimental resolution 
we convoluted the sums of weighted pDOS {'GGA + U 
DOS', solid red curves) with Gaussian functions having 
a FWHM of 100 meV and 700 meV for UPS and XPS, 
respectively (grey-shaded areas in Fig. [3]). The experi- 
mental UPS spectrum is representative of an integration 
in /c-space around the F point with a radius of roughly 
1/3 of the Brillouin zone. 
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FIG. 4: (Color online) (a) Comparison of the calculated DOS 
(grey-shaded areas) using GGA-hU with ferromagnetic spin 
alignment for different values of Uef / . While the gap between 
occupied and unoccupied states decreases from top to bottom, 
the separation between the high- and low-binding-energy re- 
gions increases. The widths of these bands, however, remain 
constant, (b) Orbital-resolved V 3d pDOS calculated using 
GGA. Three partially filled bands cross the Fermi energy lead- 
ing to a metallic state, (c) Same as in (b) but calculated with 
GGA-I-U. Only d^2_y2 and d^z axe occupied with one electron 
each, all other bands are empty. 



The basic feature observed experimentally, namely the 
separation of the spectra into a low-binding-energy part 
(ca. -2 to -4eV) of mainly V 3c? character and a high- 
binding-energy part (ca. -5 to -9 eV) derived from O 2p 
and CI 3p states is consistently reproduced in GGA-|-Ui^ 
At least for the UPS valence spectrum, the latter can 
roughly be further subdivided into a peak of mostly CI 
character for energies of approx. -4 to -6.5 eV and a struc- 
ture of mixed O and CI character below -6.5 eV. The 
changes of shape in the 0/Cl bands for different excita- 
tion energies are nicely accounted for by the appropriate 
cross-sections. While the substructure of these bands 
is reproduced remarkably well, the agreement becomes 
less satisfactory regarding the shape of the V 3d band. 
This observation will be discussed below in the context 
of ARPES data. 

An important issue when comparing PES with 
GGA-I-U (or LDA-I-U) calculations is the choice of the 
Uef f parameters for the GGA-I-U functional. In our cal- 
culations we treated only the V 3d orbitals with the Uef / 
terms. This implies that (i) the gap between the 0/Cl 
p bands and the V 3d bands closer to the Fermi edge as 
well as (ii) the charge gap between occupied V 3c? states 
and unoccupied ones will depend on the value of Ueff- 
On the other hand, the widths of the different bands are 
unaffected and remain constant upon variation of Ueff- 
In Fig. Hlja) we show the GGA-I-U pDOS for three sets 
of Ueff values. For Ueff — 1.3 eV the observed PES 
separation between 0/Cl and V occupied bands is re- 
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TABLE I: Comparison of the 3d crystal field splittings (in eV) 
from GGA calculations for VOCl and TiOCl, relative to the 
lowest 3d orbital {d^2_y2). 





dxz 


dyz 


dxy 


d,2 


VOCl 


0.025 


0.33 


1.63 


1.93 


TiOCl 


0.241 


0.46 


1.54 


2.08 



produced by the calculations (sec Fig. [3]) although for 
this value the gap between the highest occupied and the 
lowest unoccupied states is underestimated compared to 
the gap observed in photoemission as well as to optical 
absorption experiments, which reported a charge gap of 
«2eVi^i^ For f/e// = 4.3 eV the charge gap is basically 
reproduced but the 0/Cl - V gap is strongly underesti- 
mated. Since the latter is fundamentally determined by 
the hybridization pattern in VOCl, it has to be captured 
in the calculations. The parameter value Ueff = 1.3 eV is 
thus the most sensible choice for describing this system. 
This value is for instance compatible with IJ ~ 2.3 eV 
and J = 1 eV as possible onsite Coulomb and exchange 
parameters. 

The distorted nature of the VO4CI2 octahedra is re- 
flected in the crystal field splittings for the V 3(i orbital 
energies as shown in Tab. HI where these energies rela- 
tive to the lowest orbital energy are compared to those of 
TiOCl. The values were obtained by considering the first 
order moment in the DOS GGA results (see Fig. [H^b)). 
The onsite (i.e. charge neutral) orbital excitation en- 
ergies from occupied d^i^yi to the different unoccupied 
orbitals found in optical absorption experiments arc in 
good agreement with our calculated valuesi^i^i^i 

We observe that in VOCl the dj^i_yi and d^z states 
are almost degenerate unlike in TiOCl for both GGA 
and GGA-fU. Comparing the two methods, it is impor- 
tant to note that an insulating state with a finite gap at 
the Fermi energy is obtained only when correlations are 
included. This shows that correlations are the driving 
force for the suppression of metallicity in VOCl and that 
the system is a (multi-band) Mott insulator. In the case 
of TiOCl only the dxn_yi orbital is occupied, defining 
a dominant one-dimensional hopping path along the h 
directiouiii^ Nevertheless, from a quantitative compari- 
son of electronic dispersions in TiOCl and TiOBr as well 
as from an analysis of all interaction paths, the impor- 
tance of two-dimensional correlations for the electronic 
properties has been pointed out for these systems F^^i^^iii^ 

From a closer inspection of Fig. [D^c) one sees that the 
total bandwidth of dxz is larger than that of dx^-y2^ sug- 
gesting that strong antiferromagnetic exchange coupling 
would be expected along a mainly due to dxz orbitals. 
As the latter have a zig-zag overlap along the a axis 
(see Fig. [T|) one can expect that VOCl is intrinsically 
less one-dimensional from the electronic structure point 
of view compared to TiOCl, which is in line with the 
observation of a two-dimensional magnetic structure as 
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FIG. 5: (Color online) Upper panel; ARPES intensity plots 
I{k,E) of the valence region along FY (6 axis; left) and FX 
(a axis; right) together with band dispersions from GGA-f U 
{Ueff = 1.3 eV). The latter have been shifted in energy to 
account for the effects described in Sec. IIIII The similarities 
especially in the high-energy parts (O/Cl bands) are remark- 
able. Lower panel: EDCs along YFY (left) and XFX (right). 



mentioned abovc i^^'^ This is verified by calculating en- 
ergy differences within GGA+U among different mag- 
netically ordered states and then mapping the results 
to a Heisenberg model. With this method we obtain 
the exchange integrals between nearest neighbor inter- 
acting V ions along a (Ja), along b (Jf,) and along the 
diagonal direction (Jab) (see Fig. [T|). We performed the 
calculations for various Ueff values and find that Ja is 
always larger than Jb and both are antiferromagnetic 
(e.g. Ja/Jb ~ 1-1 at Ueff = 1.3 eV and Ja/Jb ~ 1-6 
at Ueff = 3.3 eV ) while Jab is smaller than Ja and 
Jb- However, the nature of the Jab exchange depends 
on the choice of the Ueff value, e.g., Jab/Jb = 0.31 at 
Ueff = 1.3eV and Jab/Jb = -0.5 at Ueff = 3.3eV. This 
is probably due to the fact that ferro- and antiferromag- 
netic contributions are equally important for the Jab ex- 
change and with Ueff increasing the antiferromagnetic 
contribution is suppressed, while the ferromagnetic con- 
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FIG. 6: (Color online) V 3d part of the valence band. Upper 
panel: ARPES intensity plots I{k, E) along FY (6 axis; left) 
and FX (a axis; right) together with band dispersions from 
GGA+U {Ueff = 1.3 eV). The latter have been shifted in 
energy to account for the effects described in Sec. IIIII Lower 
panel: EDCs along YFY (left) and XFX (right). 



tribution is barely affected, and Jab is then ferromagnetic. 
Nevertheless, the dominant exchanges J a and Jb are al- 
ways antiferromagnetic and their strengths remain com- 
parable for various C4//, indicating that VOCl should 
be a two-dimensional antiferromagnet at low tempera- 
tures without a preferred onc-dimcnsional hopping path 
induced by largely different hopping matrix elements as 
in TiOCl. Indeed, the NMTO downfolding results iden- 
tify the dominant hopping parameters to be of the same 
order of magnitude in a (0.18 eV) and h (0.13 eV) di- 
rections and a factor ~ 0.6 smaller in the c direction. 
This qualitative trend in dominant hoppings agrees well 
with the estimation of the exchange parameters as given 
above. 

In order to elaborate some more on the dimensionality 
issue we compare ARPES data with band structures from 
GGA-hU in Figs. and [51 First, Fig. [5] shows intensity 
plots of the dispersions of the entire valence band in the 
Brillouin zone along FX and FY corresponding to the 
crystallographic a and b axis, respectively, together with 
the respective energy distribution curves (EDCs). Also 
shown is the band structure calculated using GGA-I-U. 

A large number of different bands and rich dispersive 
features which differ between the two directions arc ob- 



served for the 0/Cl bands in very good agreement be- 
tween theory and experiment. This also indicates the 
good quality of our single crystals. In the low-binding- 
energy region, GGA-fU shows a total of four V 2>d bands 
since there are two electrons per V site and two V sites 
per unit cell. In the ARPES data, however, the V id 
weight appears as a rather broad and weakly dispers- 
ing hump without any sharp features, as can be seen 
in the blow-up shown in Fig. [6l This behavior was re- 
producibly observed on many different samples and is 
not due to, e.g., a limited experimental resolution. It 
is rather an intrinsic property which for several reasons 
cannot be expected to be fully captured in our calcula- 
tions. From a strict many-body view the V 3c? bands 
do not have quasi-particlc character, but must be in- 
terpreted as incoherent weight in the spirit of a Mott- 
Hubbard scenario. Thus, a description of the detailed line 
shapes within GGA-I-U, being an effective single-particle 
approach which does not account for dynamical and spa- 
tial fluctuations, is admittedly intricate. This issue has 
been extensively studied for the case of TiOCl, where 
dynamical mean field theory (DMFT) with a Quantum 
Monte Carlo (QMC) impurity solver— in order to take 
full account of dynamical correlations, as well as the ex- 
tension to cluster- rather than single-site calculations;^ 
have led to considerable improvement. Of course, due to 
the predominant one-dimensionality in TiOCl as well as 
its strong tendency to form singlets (evidenced by the un- 
conventional spin-dimerizcd phases at intermediate and 
low temperatures^) the influence of such fluctuations is 
expected to be stronger than in VOCl. Nevertheless, this 
issue calls for further investigation in the present case as 
well. 

Despite these qualifying remarks, we shall give a more 
detailed analysis in order to compare to the computa- 
tionally equally complex case of TiOCl. Compared to the 
latter, the differences of the ARPES dispersions of the V 
3d weight between the a and b axis are less pronounced^ 
with the peak maximum appearing at highest binding en- 
ergies at the F point and shifting slightly to lower bind- 
ing energies towards the zone boundary in both cases. A 
comparison shows that the trend towards lower binding 
energies as well as the evolution of the spectral weight dis- 
tribution are qualitatively in line with the calculations: 
large weight is found along a half way between F and X, 
along b right at the Y point, i.e., where bands cross in 
both cases. The shape of the peak, which can be seen in 
the lower panels of Fig. [SI reflects these trends as well: 
along b it develops from a flat hump to a more distinct, 
symmetric peak towards the zone boundary, while along 
a the slight asymmetry towards the high-binding-energy 
side at F shifts towards the low-binding-energy side at 
X. The overall width at the X point appears to be larger 
than at the Y point, in line with the calculated result 
that the four GGA-I-U bands merge at Y while two pairs 
remain well separated at X. However, a detailed inter- 
pretation of this behavior needs further analysis which is 
beyond the scope of the present work, although we want 
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to note that the observed asymmetric shape in VOCl is 
reminiscent of TiOCl where the asymmetry has been at- 
tributed to correlations!^^ 

Overall, the agreement between GGA+U and ARPES 
for VOCl is better than in the case of TiOCl. This is 
probably due to the fact that the influence of dynamical 
and/or spatial fluctuation effects in VOCl is smaller than 
in TiOCl. 



V. CONCLUSIONS 

In conclusion, we have investigated the electronic 
structure of the low-dimensional 3(f system VOCl using 
photoemission spectroscopy and GGA-hU band structure 
calculations. From a comparison to GGA calculations 
we could show that its insulating nature is due to strong 
correlations of the V 3d valence electrons instead of band 
structure effects, which is the first conclusive evidence 
for the Mott character of the insulating state. The issue 
of choosing Uef / = U — J has been addressed concern- 
ing the correct description of the hybridization pattern 
in favor of a reasonable magnitude for the charge gap. 



Calculation of the exchange integrals shows that VOCl 
is two-dimensional both electronically and magnetically, 
and the obtained antiferromagnetic exchange within the 
ab plane is robust upon variation of Ueff and in line 
with experimental observations. Calculated band disper- 
sions agree well with the ones observed via angle-resolved 
photoemission spectroscopy for the high-binding-energy 
uncorrelated bands, and qualitatively match the V 3d de- 
rived weight. No considerable in-plane anisotropy within 
the VO bilayers is observed in ARPES, supporting the 
two-dimensional nature of the system. Overall, we could 
show that VOCl is a quasi-two-dimensional (as opposed 
to quasi-one-dimensional) system and can be well de- 
scribed as a two-band Mott-Hubbard insulator with d^ 
ground state configuration. 
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